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The  coefficient  of  compatibility  between  two  menu  items  is  defined  as  the 
difference  between  the  conditional  and  marginal  probabilities  of  selecting  one  after 
the  other.  Maximum  likelihood  estimates  from  selection  statistics  are  given,  and 
the  asymptotic  distribution  of  the  estimator  is  determined  so  that  test  of  null 
hypothesis  and  confidence  intervals  regarding  the  coefficients  can  be  obtained.  A 
method  is  described  for  building  data  base  for  large  joint  probability  matrices  from 
menu  selection  data. 
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STATISTICAL  ESTIMATION  OF  THE  COEFFICIENT 
OF  PAIRWISE  COMPATIBILITY  FOR  MENU  ITEMS 

by 

Joseph  L.  Balintfy 
and 

Ram  C.  Dahiya 

I  INTRODUCTION 

With  recent  advancements  in  the  art  of  representing  food  and  meal 
preferences  by  mathematical  models  [  2  J  along  with  the  possibilities  of 
defining  and  solving  optimum  human  diet  problems  by  mathematical  program¬ 
ming  techniques  [  1  ],  difficult  questions  emerge  concerning  the  role  and 
measurability  of  compatibility  between  menu  items. 

A  meal  selected  by,  or  planned  for  an  Individual  usually  consists  of 
a  set  of  menu  items  which  are  mutually  complementary  in  the  sense  that  each 
item  represents  one  of  the  courses  of  the  meal.  It  is  tempting  to  consider 
the  utility  of  a  meal  in  terms  of  the  utilities  of  its  components.  Recent 
results  in  multiattribute  utility  theory  open  the  way  for  a  variety  of 
representations,  and  this  is  where  the  issue  of  compatibility  comes  in. 

The’ simplest  additive  utility  model  would  imply  that  the  utility  of  a 
meal  is  equal  to  some  weighted  sum  of  the  univariate  utilities  of  the  menu 
items  in  the  meal.  Additivity,  therefore,  means  that  the  utility  of  the 
meal  is  completely  explained  by  the  utilities  of  the  items  irrespective 


of  their  relative  combination.  Empirical  evidence  I  5  ]  is,  however,  over¬ 
whelming  that  this  is  not  so.  People  find  some  combinations  of  items  more 
or  lees  compatible  than  others,  meaning  that  compatibility  is  also  a  factor 
in  the  utility  of  a  meal. 

If  the  set  of  menu  items  under  consideration  is  preference  and  there¬ 
fore  utility  independent,  a  measure  of  comparability  can  be  derived  from 
the  multiattribute  utility  model  of  Keeney  [6  ].  This  model  expresses  the 
utility  of  a  meal  as  a  weighted  sum  of  additive  utilities  plus  a  weighted 
sum  of  all  possible  crossproducts  of  univariate  utilities.  Depending  upon 
the  sign  of  the  coefficients,  this  second  sum  may  increase  or  decrease  the 
value  of  additive  utilities,  and  hence  would  represent  the  effect  of  compat¬ 
ibility  in  a  given  combination  of  items.  Unfortunately,  again,  utility 
independence  can  not  be  assumed  to  be  true  for  any  set  of  menu  items,  and 
not  even  for  any  set  of  pairwise  combinations  of  items.  Consequently,  re¬ 
search  is  still  in  progress  to  find  the  appropriate  expression  of  compati¬ 
bility  in  utility  terms. 

In  a  recent  report  by  Balintfy  and  Sinha  (  3  }  a  meal  planning  model  is 
presented  which  does  not  require  the  representation  of  compatibility  of 
menu  items  in  the  utility  measure  of  the  meals.  The  idea  is  advanced  that 
a  mathematical  programming  model  which  maximizes  an  additive  utility  function 
could  determine  which  items  would  be  most  liked  by  a  given  population  in  a 
given  time  period,  and  from  this  fixed  set  of  items  a  separate  scheduling 
activity  could  combine  the  items  into  a  sequence  of  meals  on  the  criterion 
of  compatibility  alone.  This  way  the  concepts  of  the  utility  and  compati¬ 
bility  of  a  meal  are  conveniently  separated;  the  first  being  included  only 


in  the  objective  function  additively,  while  the  second  appears  among  the 
scheduling  constraints,  and  hence  is  no  longer  linked  directly  to  utility 
measures.  Consequently,  it  can  be  represented  by  any  measure  for  which  data 
collection  is  feasible. 

This  paper  suggests  the  introduction  of  a  statistical  measure  for 
compatibility  of  menu  items  in  the  form  of  a  coefficient  of  pairwise  compat¬ 
ibility  of  menu  items.  This  coefficient  is  defined  on  a  cartesian  product 
set  of  two  sets  of  menu  items  where  the  sets  under  consideration  are  ident¬ 
ified  with  two  nonidentical  courses  of  the  meal,  such  as  entrees  and  veget¬ 
ables.  Limiting  the  notion  of  compatibility  to  pairwise  compatibility  alone 
is  arbitrary,  but  it  is  the  first  step  of  the  investigation  which  may  lead 
to  further  extensions  later. 

The  first  part  of  the  paper  provides  the  definition,  some  properties 
and  examples  of  the  coefficient  of  pairwise  compatibility.  The  second  part 
describes  the  method  of  estimation  and  the  probability  distribution  of  the 
estimator.  The  last  part  is  devoted  to  the  techniques  of  estimating 
coefficients  for  the  whole  cartesian  product  set  from  subsets  of  data  as 
it  is  usually  available  through  a  sequence  of  selective  menu  schedules. 

The  application  of  this  coefficient  to  menu  scheduling  algorithms  and  ex¬ 
perience  with  data  collection  will  be  the  subject  of  later  reports. 


IX  PROBLEM  DEFINITION 


Let  X  represent  a  finite  set  of  menu  items  such  as  entrees,  and  let 
Y  represent  another  finite  set  of  items  such  as  vegetables.  Then  the 
Cartesian  product  set  U  =•■  [X  x  Y]  represents  all  pairwise  combinations  of 
menu  items  under  consideration.  The  elements  of  U  can  be  considered  as 
elementary  meals  consisting  of  only  two  courses.  If  x^  is  one  of  the 
entrees,  i.e.,  X|  e  X  ,  and  yj  is  the  j-th  vegetable,  i.e.,  y^  e  Y  , 
then  the  pairwise  combination  of  items  (x^,yj)  denoted  by  u-jj  is  the 
element  of  U  ,  i.e.,  (x-^y-j)  ■  uy  b  U  . 

If  one  considers  sets  X  and  Y  as  choices  of  menu  items,  such  as 
entrees  and  vegetables  respectively,  then  it  is  conceivable  that  an  indi¬ 
vidual  will  have  a  linear  preference  ordering  at  any  point  of  time  for  all 
combinations  (x^,yj)  which  is  based  on  his  utility  for  x^  ,  his  utility 
for  yj  and  his  utility  for  the  compatibility  between  and  y_j  in  a 

meal.  For  a  rational  person,  this  preference  ordering  is  revealed  by  his 
choice-behavior  in  selecting  a  particular  combination  of  items  in  a  sequence 
of  trials.  Consequently,  one  can  assign  probability  measures  for  the  sets 
X  ,  Y  and  U  as  follows: 


P(Xi)  >.  o  j 

2  pCx*)  *  1 
x^eX 

’  p(yj)  1  o  ; 

2  p  Cy.* )  -  1 
yj£Y 

p(xi»yj)  >  o 

;  2  pCx^y-j)  *  1 

uijeU 
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For  the  sake  of  notational  simplification  p(xi)  =  pi  ,  p(y.)  =  p.  and 

P (xi>yj )  =  Pij  notations  also  will  be  used  in  subsequent  parts  of  the  paper. 
The  first  two  of  these  notations  will  be  referred  to  as  the  marginal  proba¬ 
bility  of  selection,  while  the  last  term  is  called  the  joint  probability  of 
selection.  These  terms  are  related  by  the  well  known  identities 

p(xx)  =  Z  pCx^y.) 
yj£Y 

p(yj)  *  £  p(xi>yj) 

xieX 


Decisions  in  choosing  any  item  combination  from  set  U  are  usually  made 
sequentially  in  the  sense  that  a  choice  is  made,  say  from  set  X  first 
(entrees) ,  followed  by  a  choice  from  set  Y  ,  implying  that  the  choice  of 
yj  may  be  conditioned  by  the  choice  of  a  particular  x±  .  Consequently, 
one  can  talk  about  the  conditional  probability  of  selecting  y.  given  that 

•J 

x^  is  already  selected  according  to  the  formula 


pCyj  Ixi) 


p0Vyj) 

PCXi) 


which  means  that  the  conditional  probability  of  selections  is  completely 
defined  by  the  corresponding  marginal  and  joint  probability  measures.  Fur¬ 
thermore,  it  is  also  known  that  if  pCx^y.)  =  p(x^)p(y.)  ,  then 

J 

P(Yj  Ixf)  -.=  P(Yj)  >  which  means  that  the  selection  of  y^  is  unconditional 
of  the  choice  of  x^  . 

It  is  proposed  that  a  coefficient  of  pairwise  compatibility  of  menu 
items  be  created  on  the  basis  of  the  above  discussed  conditional  dependence 


-5- 


criterion.  Accordingly  a  coefficient  is  defined  as 

(1)  ©ij  =  p(yj  \x±)  -  p(yj) 

to  express  the  compatibility  between  menu  items  x.^  and  y^  according  to 
the  following  rules  of  interpretation: 

item  x^  is  compatible  with  item  y^  , 
item  x^  is  not  compatible  with  item  y^  . 
the  items  are  neither  compatible  nor  incompatible,  hence 
they  are  utility  additive. 

It  is  noteworthy  to  add  that  the  above  definition  of  the  0^  coefficients 
can  express  degrees  of  compatibilities  or  the  lack  of  it  as  well  as  statis¬ 
tically  significant  differences  from  zero. 

Example:  Consider  a  set  of  three  entrees,  X  =  {x  ,x  ,x  }  and  three  veg- 
etables,  Y  =  (y^y 2,y 3}  .  Let  P  =  |  |p-^j|  |  be  the  matrix  of  joint  proba¬ 
bilities  and  p (x)  and  p(y)  the  vectors  of  marginal  probabilities  over  X 
and  Y  respectively.  Let  some  arbitrary  values  of  these  probabilities  be 
as  follows: 


i/7 

11/3? 

1/18 

1/6 

1/9 

p(x)  = 

1/3 

.  ;  p  Cy)  = 

14/36 

;  p  = 

5/36 

1/9 

1/12 

_}(  3_ 

11/36 

JV9 

1/9 

i/9_ 

where  the  elements  p^  indicate  the  joint  probability  in  the  i-th  row  and 
the  j-th  column. 

From  these  data  by  formula  (1)  the  matrix  of  conditional  probabilities 


0ij  > 0 
®lj  <  0 

©ij  »  0 
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PCyj lxi)  denoted  by  Py|x  can  be  computed,  and  yields 


Then  subtracting  p(y)  from  each  column  of  P  i  yields  the  matrix 

y|x 

of  the  coefficients  of  compatibility,  ll^ijll 


-5/36 

4/36 

1/36 

0ijll 

= 

4/36 

-2/36 

-2/36 

1/36 

-2/36 

1/36 

It  is  interesting  to  note  that  the  coefficient  0^  will  take  negative 

and  positive  values  even  in  the  last  row  of  the  table  where  the  input  data 

would  imply  indifference.  One  reason  for  this  is  that  pO^yj)  t  p(x3)p(y.). 

Furthermore,  just  because  (x^),  (x^),  (x^)  are  equiprobable  does  not 

mean  that  (x^)  are  indifferent.  In  fact,  if  one  looks  at  the  choices  of 

Xl  and  x^  combined,  then  one  finds  that  y^,y3  are  equally  preferred, 

but  y2  is  preferred  over  both  y^  and  y3  .  Now  just  the  fact  that  y3 

gets  equal  weight  to  y^  and  y2  in  combinations  with  x3  means  that  y3 

does  not  go  well  with  xg  and  this  is  indicated  by  the  number  -2/36.  The 
■» 

above  paradox  can  be  generalized  as  follows.  Suppose  p(ya),  p(y2),  •**, 

P(yj)  are  the  marginal  probabilities  of  y1 ,y2 , • • • ,yj  found  from  (1-1) 
entrees.  Now  suppose  the  I-th  entree  is  selected,  so  that  pCy^Xj)  =  1/J 
for  all  j  uniformly.  Also,  let  max  pCy^)  =  p(ym)  ^  1/J.  Then  this  implies 
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that  the  m-th  vegetable  gives  a  worse  combination  with  x-j-  .  Also?  if 
min  P(y-j)  =  p(Yn)  >  then  this  implies  that  yn  goes  best  with  the  I-th 
entree.  Hence  we  will  still  get  negative  and  positive  values  for  indicating 
compatibility  of  the  I-th  entree  with  different  vegetables. 

The  above  example  also  shows  that  the  coefficients  of  compatibility 
depend  intimately  upon  the  cardinality  of  the  sets  of  items  involved  in 
defining  the  measure.  Therefore,  for  the  purpose  of  finding  coefficients 
which  are  invariant  under  various  conditions  of  application,  the  cardinality 
of  U  should  be  as  large  as  possible.  Since  choices  satisfying  this  con¬ 
dition  are  rarely,  if  ever,  available  and  hence  observable  on  selective 
menus,  a  procedure  is  outlined  in  the  last  section  of  the  paper  to  allow 
the  construction  of  the  joint  probability  matrix  P  from  estimates  of  its 
submatrices. 
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Ill  THE  ESTIMATION  OF  THE  COEFFICIENTS  OF  COMPATIBILITY 


In  order  to  estimate  0^^  ,  defined  in  section  II,  we  start  with  n 
independent  meal  selections.  If  we  want  to  estimate  0^  for  a  population, 
then  n  is  a  random  sample  of  people  belonging  to  this  population.  On  the 
other  hand,  n  could  be  different  meal  selections  by  the  same  individual 
if  one  is  interested  in  estimating  for  an  individual.  Let 

rij  =  number  of  meal  selections  (out  of  n)  containing  (x^,yj). 

ri  =  number  of  meal  selections  (out  of  n)  containing  menu  item  x^  . 

rj  =  num^er  °f  meal  selections  (out  of  n)  containing  menu  item  y-  . 

J 

i  =  1,  2,  •••,  I  ;  j  =  1,  2,  ...,  J  . 

Then  the  maximum  likelihood  estimator  of  0^  is  given  by 

(2)  -  rj/n 

The  estimator  0^  is  an  unbiased  estimator  of  0^  ,  since 


where  E  denotes  expectation  with  respect  to  r .  .  Derivation  of  (3  ) 
ri  , 

uses  the  fact  that  r^j |r^  is  a  binomial  random  variable  with  parameters 


Now  we  develop  the  asymptotic  distribution  of  0^  which  is  required 
for  testing  the  hypothesis  that  O^j  =  0  .  This  distribution  can  also  be 
used  in  order  to  find  the  confidence  interval  of  0^..  , 

Theorem 

The  asymptotic  distribution  of  ^  (0^-0^)  is  N(0,a^)  ,  where 


fij  ~  pij(1  "  Pij/Pi  +  2Pij  "  2Pi)/Pi  +  Pj(l~Pj) 


Proof 


^(§ij 


-  0ij)  =  /n  { 


=  ^  rfii  _  D  _  ZiL 


n  Pi  + 


n 


pj) 


where 


°n  Un  “  vn 


n 

ri 


UT 


VT 


&  (hi  _  p±i  ) 

n  n  Pl 

^  i1  ~  pj  > 


Since  (Un,  Vn)  have  a  limiting  distribution  and  aR  - >  ^  ,  making 

use  of  the  result  in  Rao[  7>p.319],  the  asymptotic  distribution  of 
/n  (0ij  "“©ij)  is  the  same  as  the  asymptotic  distribution  of  Tn  where 

<*>  Tn  ■  FiUn-Vn 

Now  in  order  to  obtain  the  asymptotic  distribution  of  Tn  ,  we  define 
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Then 


1  if  k-th  meal  selection  contains  (x.,y.) 
Xijk  =  1  J 

0  otherwise 

lc  =  1 ,  2 ,  •  •  • ,  n 


T  =  y  /v  fii  v 

n  Pi  kfi(Xi^k  “  Pi  Xi*k  "  Pix-jk  +  PiPj )/n 


where 


_  n 

=  A  Z  Yijk/n  , 
k=l 

J 

Xi‘k  =  xijk  > 


X-jk 


2  xijk  , 
i=l  J 


Yijlc  -  (xijk  "  Xi*k  -  Pi  X.jk  +  piPj)/Pi  . 


We  need  the  first  two  moments  of  Yijk  in  order  to  obtain  the  asymP 
totic  distribution  of  Tn  .  It  is  straightforward  to  show  that 

E(Tn)  =  0 

In  order  to  obtain  V(Yi:jk)  ,  we  need  the  following  results. 

E(xijk)  =  Pij.  , 

** 

V(Xijk)  =  Pijqij  >  where  q...  =  1-Pij  , 
c°v (xijk>  X^jk)  =  -PijP£j  i  +  Z  , 

V(X.jk)  =  Pjqj  ,  q.  =  1-Pj  , 


-11- 


vcxt.k)  =  Pim  , 


Cov(XljkX.jk) 


Cov  (Xijk,x1-k) 

Cov(Xi.k,X.jk) 


I 

X  Cov(XljkXjl.k) 

1 

I 

PijCl-Pij)  -  Z  PijPAj 
&—  1 
l+x 

Pijqj  > 

Pijqi  > 

E(xi.kX.jk)  -  ECXi.k)  E(X. jk) 
p(xijk=1)  -  PiPj  > 


=  Pij  -  PiPj  • 


Making  use  of  the  above  results  in  the  following  equation, 

Pi  v<Yijk)  =  v^xijk)  +  v(xi-k>  +  Pi  V(X.jk)  -  2  CovCX^^.k)  - 
"  2  Pt  c°v (xijk>x* jk)  +  2  Pij  Cov(X^.k,x. jk>  > 

and  simplifying,  we  have 

(5)  a±j  =  V(Y-^jk)  =  PijCl-pij/pi  +  2Pij  -  2Pi)/p?  +  pjqj 

Since  Y-jjk  ,  k  =  1,  2,  ••*,  n  ,  are  independent  and  identically 
distributed  random  variables  with  mean  zero  and  variance  a.,  given  by  (5  ), 

•J 

by  Central  Limit  Theorem  and  equation  ( 4  )  we  have 

d  2 

Tn  - >  N(0 ,ai:j)  , 


-12- 


where  - >  denotes  convergence  in  distribution.  Consequently,  we  have 

proved  the  result  in  Theorem  that 


-  By)  - >  N(0,a?-j)  . 


Now  in  order  to  obtain  confidence  internal  of  0^  ,  we  have  the 
following  Lemma  which  can  easily  be  proved  by  making  use  of  result  in 
Cramer  [  4  ,p254] . 

Lemma  1 : 


The  asymptotic  distribution  of 


A(°ij  "  0ii) 
\  3 


is  N(0,1)  ,  where  a 


ij 


is  an  estimator  of  obtained  on  replacing  by  r^/n  ,  p-^  by 

r^/n  and  pj  by  rj/n  in  . 

The  test  statistic 

Hj  "  ^ 

can  now  be  used  for  testing  the  null  hypothesis 


Ho  :  ©ij  "  0  , 

where  2^.  has  a  standard  normal  distribution  for  large  n  . 
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IV  OBTAINING  JOINT  PROBABILITY  ESTIMATES  FROM  SUBSETS  OF  THE  FULL  PRODUCT  SET 


It  is  obvious  that  the  estimation  of  the  coefficient  of  compatibility, 

©ij  »  any  pair  of  menu  items  is  related  to  the  estimation  problem  of  the 
corresponding  joint  probability  term  p^  which  in  turn  is  the  function  of 
the  cardinality  of  the  product  set  [X  x  Y]  .  Ideally,  therefore,  sets  X 
(entrees)  and  Y  (vegetables)  should  be  large  enough  to  include  all  the  items 
an  individual  may  routinely  encounter  in  his  menu  selection  activity.  The 
problem  here  is,  of  course,  that  the  number  of  choices  which  can  be  offered 
at  any  time  in  a  selective  menu  are  much  less  than  the  elements  in  X  or  Y  . 
Consequently,  observations  are  feasible  only  on  subsets  of  [X  x  Y]  and  a 
technique  is  necessary  to  project  the  joint  probabilities  of  the  subsets  into 
those  of  the  full  set  under  consideration.  We  will  show  first  that  the  reverse 
process  is  always  feasible,  and  from  that  the  procedure  of  correct  projection 
easily  follows. 

Let  matrix  P  =  |  Ip^ | |  denote  the  joint  probabilities  of  selection 
when  sets  X  (entrees)  and  Y  (vegetables)  are  of  cardinality  I  and  J 
respectively,  and  when  the  cartesian  product  set  U  =  (X  x  Y)  represents  the 
full  set  of  pairwise  combinations  of  all  the  items  under  consideration.  Let 
Sk  be  a  subset  of  U  ,  and  consider  K  such  subsets  where  U  for 

k  =  1,  2,  •••,  K  holds.  The  subset  itself  is  the  cartesian  product 

set  of  item  sets  X^  and  Yjc  with  cardinality  1^  <  I  and  <  J  respect¬ 
ively.  For  each  subset  Sk  there  exists  a  corresponding  submatrix  of  the 
joint  probabilities  P  denoted  by  P^. 

Now  consider  the  joint  probability  matrix  P(Sk)  defined  only  on  the 
limited  choices  X^  and  Y^  .  Clearly,  P^  ^  P(Sk)  ,  but  they  are  related 
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by  a  Simple  rule  of  proportionality  If  one  makes  the  llgltlmate  assumption 

that  the  preference  ordering  of  the  individual  remains  unchanged.  In  this 
case 

(6)  P(Sk)  =  Pk/ak 
where 

Ik  Jlc 

ak  =  *  *  Pij 

i=l  J=1 

and  is  the  element  of  the  known  matrix  Pk  .  The  constant  ak  is 

nothing  else  but  the  sum  of  the  probabilities  in  matrix  Pk  which  is  also 
equal  to  the  sum  of  the  marginal  probabilities  for  the  given  submatrix  with 
respect  to  subsets  Xk  or  Yk  and  hence  ak  is  less  than  one. 

Now  it  is  obvious  that  we  can  project  PR  into  P(Sk)  ,  but  our  problem 
is  to  go  backward,  because  in  reality  only  estimates  of  P(Sk)  will  be  avail¬ 
able,  and  the  Pk  submatrices,  and  eventually  the  matrix  P  will  have  to  be 

K 

u  =  (J  Sk  in  such  a  way  that  the  K 
k=l 

K-l  nonempty  intersections  of  the 
such  an  intersection  set. 


reconstructed  by  some  process. 

First  we  make  the  assumption  that 
subsets  are  not  disjoint,  but  at  least 
type  Skfl  S£  exist.  Let  ^k£  denote 
Then  we  have  from  (6) 

z  Pij(sk)  =  2  Pij/ak 

’  kA  °k* 

and 

1  pij(V  =  1  Pij/af 

V  °k* 
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Consequently 

Va*  -  z  Py(s*)/  z  Py(sk)  =  bk, 

°k*  “id 

or 

ak  ~  bk£a£  =  0 

where  b^  is  considered  known,  and  afc  and  are  to  be  determined. 

If  we  consider  the  K-l  intersections  between  the  K  subsets,  and 
make  the  k+1  =  *  substitution,  we  obtain  a  system  of  K-l  equations 

ak  “  ckalc+l  =  0  ;  k  =  1 ,  2 ,  •  •  • ,  K-l 

where 

ck  =  bk,lc+l 
Also,  we  have 

ax  +  a2  +  • • •  +  aR  =  C 

This  system  has  a  unique  solution  for  any  given  value  of  C  .  In  matrix 
vector  form  we  have 


which  can  be  stated  equivalently  by  the  matrix  algebraic  symbols 
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(7)  Ma  =  d 

where 

a  and  d  are  column  vector  symbols  for  the  variable  a,  •  k  =  1  ?  . . . 

and  for  the  constant  vector  at  the  right  hand  side.  The  solution  of  (7)  can 
be  written  as 

a  =  H_1d 

which  is  equivalent  to 

a  =  cmk 


where  M}c  denotes  the  K-th  Clast)  column  of  M_1.  Because  of  the  special 
structure  of  M  the  components  of  ^  ,  and  hence  that  of  a  ,  can  be 
explicitely  computed  as  follows: 

,/A 


=  C.c  .c  . 
1  12 


'K-1‘ 


where 


a2  =  C.c2.c3...-cK_1/A 


lK-l  =  C-ck-1/A 


aK  =  C/A 


in  this  case. 
Since 


A  =  |M|  = 


1  +  CK-1  +  CK-2 


CK-1  + 


+  c  .c 
1  2 


*  *  C 


K-l 


a  unique  solution  can  be  found  for  any  value  of  C  .  The  correct  value  of 
C  obtains  by  normalizing  the  resulting  P  matrix  elements  to  satisfy  the 
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1 


I  J 

z  z 

i=i  j=: 

conditions  for  probability  measures. 

The  above  derivations  open  the  way  for  experiments  and  experimental 
designs  with  the  objective  of  obtaining  estimates  of  P  from  observations 
available  on  a  sequence  of  selection  statistics  from  selective  menus.  It 
is  believed  that  the  procedures  outlined  are  necessary  and  sufficient  for 
the  statistical  basis  of  establishing  food  compatibility  data  banks  for 
computerized  food  service  organizations. 
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